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Abstract 

We propose a scenario where the pattern of the recently observed charge-dependent az- 
CN imuthal correlations can be understood qualitatively. This is based on the cluster picture and 

O the assumption that the charged hadrons that flow outward from the surface of the overlapping 

^ region of the colliding nuclei move primarily parallel to the reaction plane. We also point out 

the there is a strong electric field induced by the transient magnetic field during the parton 
production in the initial phase of the relativistic heavy-ion collision and discuss its possible 
relevance to the scenario. 
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cn 1 Introduction 

oo 

oo 

Relativistic heavy-ion collision provides a fertile ground for a subtle interplay between electrody- 
0^ namics and QCD, since the magnetic field is very strong and short-lived so that a strong electric 
2 field is induced according to Faraday's law. Recently, the charge-dependent azimuthal correlations 
T-H were detected by STAR at Relativistic Heavy-Ion Collider p[] in an attempt to detect the chiral 
>. magnetic effect (CME) due to a possible local strong CP violation [21 |3j. Since all the direct 
*^ parity-odd observables would vanish due to the fact that the induced current in CME can either 
S-H align or anti-align with the angular momentum, the STAR experiment measured the correlator 
(cos(0Q, + 4>/3 — 2\l//jp)), where (pa and (pp are the azimuthal angles of the charged particles and "^rp 
is that of the reaction plane. This quantity is not parity- violating [H, H], but can be sensitive to 
parity- violating effect. It is shown by STAR that in the three-particle correlator, the explicit 
determination of the reaction plane is not needed [5] . Assuming that particle c in the three-particle 
correlation is correlated to particles a and only in common correlation to the reaction plane, one 
has (cos(0Q, + 4>i3 — 20c)) = (cos(0a + 0^ — 2"^ Rp))v2,c, where V2,c is the elliptic flow of the particle c. 

Based on this assumption, the experimental three-particle azimuthal correlation and the sepa- 
rately measured correlation (cos(0o — 0/?)) are combined to obtain the correlations of products of 
the azimuthal angles (sin(0a) sin(0^)) and (cos(0q,) cos(0^)) for the like-sign charged pairs (+,+), 



(-,-) and unlike-sign charged pairs (+,-)) (->+) [6]- The frame is chosen in such a way that the 
reaction plane is in the x — z plane and the y axis is perpendicular to the reaction plane. We shall 
adopt the same frame which is illustrated in Fig. [l] In this case, both the angular momentum L 
and the central magnetic field B are along the the positive y axis. The salient features revealed 
from this analysis |6] for these charge-dependent azimuthal correlations are the following: 

• For the like-sign (LS) charged pairs, it is found 

{sm{(j)^)sm{(j)fi))Ls -0, (1) 
(COS(0„) COs{4>fj))LS < 0. (2) 

• For the unlike-sign (US) charged pairs, on the other hand, one finds the following relation 

(sin(0„) sin(0^))(;5 ~ (cos(0o) cos{(f)fs))us > 0. (3) 

Furthermore, the smallness of (sin(0Q,) sin(0/3))2,s' in Eq. ([T]) holds for the range of centrality from 
65% (peripheral collision) down to ~ 0% (cental collision); while (cos(0a) cos(0;3))l5 is increasingly 
more negative from central to peripheral collision. On the other hand, (sin(0Q,) sin(0^))[75' and 
(cos(0Q,) cos(0/3))^5 for the unlike-sign cases are increasingly more positive from central to peripheral 
collision. They are almost identical from 0% to ~ 60% centrality and are about a factor of two 
larger in magnitude as compared to (cos(0q,) cos(0/3))l5' (Eq. (|2])) in this centrality range. 

The smallness of (sin(0Q,) sin(0^))i5 is taken to imply that the like-sign charged pairs are 
mostly in plane, since sin(0) is non-negative for all the azimuthal angles [B]. The fact that 
(cos(<^Q,) cos((/)/3))l5 is negative is taken to suggest [6] that the the like-sign pairs are not only 
in plane but more back-to-back than in the same x-direction (N.B. the azimuthal angle is mea- 



sured from the reaction plane x — z in reference to positive x (See Fig. 2(a)). This seems to be 
in contradiction with CME which predicts that the like-sign pairs are out of plane. However, the 
signals of CME could be smaller than the background events due to the strong electric and mag- 
netic fields and the omnipresent strong interaction. In this case, it is essential to understand the 
background so that one can hope to pull out the signal of CME. In the case of unlike-sign pairs, 
the positive (cos(0a) cos(0^))(75 suggests that they are more in the same x direction, i.e. with the 
relative azimuthal angle 1 0^ — 0^1 < 90°. However, it is hard to comprehend why (sin(0Q,) sin(0^));75' 
and (cos(0a) cos(0/3))[/5' should be the same [H]. Since there is an elliptic flow, one would expect 
(cos(0Q,) cos(0/3))c/5 which favors the in-plane flow to be larger than (sin(0Q,) sin(0/3))(75' which fa- 
vors out-of-plane flow. Moreover, it is not clear as to why the sizes of (sin(0Q,) sin(0/3))(75' and 
(cos(0Q,) cos{(f)i3))us are a factor of two larger than that of (cos(0q,) cos(0^))l5, with the latter being 
negative. 

We shall present a scenario which may allow us to understand the azimuthal correlations of 
charged hadrons in a qualitative manner. The observation of clustering and its implications will be 
discussed in Sec. |2} The time dependence of the magnetic field and the induced electric field in the 
initial phase for the nuclear collision and their impacts on the initial parton momentum direction 
are presented in Sec. |3j Partly motivated by the large electric field present in the early phase of the 
collision, we postulate in Sec. |4] that the azimuthal correlations mainly originate from the clusters 
on the surface of the overlapping region and proceed to calculate the charge-dependent azimuthal 
correlations. We end with a conclusion in Sec. [H 
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2 Clustering 



We first address the observation that the azimuthal correlations in Eqs. ([2]) and ([3]) increase from 
the more central colhsions to peripheral ones, i.e. as % of centrality increases. This can be largely 
understood as due to clustering. In measuring the rapidity and azimuthal correlations of two charged 
particles produced in the pp collision in the central rapidity region, it was observed [S] that if one 
multiplies n — 1, where n is the number of charged particles, to the azimuthal correlations, they 
turn out to be independent of n within errors. This is an indication of short-range correlation which 
can be described by independent emission of low-multiplicity clusters. The n — 1 dependence comes 
from the fact that the correlation scales with n which is proportional to the number of clusters, 
while the total charged particles scales with n{n — 1). The same analysis based on the cluster model 
was carried out for Cu + Cu and Au + Au collisions at center of mass energy of 200 GeV p^ . 
On top of the n — 1 scaling observed in the pp case, there is an additional suppression of effective 
size (multiplicity) of the cluster toward more central collision, and more so with the away-side 
(90° < A0 < 180°) pairs than with the near-side (0° < A0 < 90°) pairs. It is suggested that 
in more central collisions, the correlated decays from the clusters, in space-time, inside the dense 
overlapping region get lost and only the correlations from the clusters on the surface survive. Thus, 
it is the combined effects of n — 1 scaling and the survival of correlations from the clusters on the 
surface which are responsible for the observed decrease of the azimuthal correlations toward more 
central collisions, as the centrality is measured by the inelastic cross section or the total multiplicity 
of particle production. 

Following this natural explanation of the centrality dependence, we shall explore the possibility 
of understanding the azimuthal correlations in Eqs. ([l])-(|3| in terms of the cluster description where 
the correlated pairs emerge from within those clusters which are mostly on the surface of the 
overlapping region. The effects of cluster particle correlations on two- and three-particle azimuthal 
correlations and their parameters have been studied from the STAR experimental input [12]. We 
shall present a scenario based on the cluster model to explain the azimuthal correlations themselves. 



3 Electric and Magnetic fields 

Next, we consider the time dependence of the electric and magnetic field landscape in the initial 
phase of the nuclear collision at relativistic speed as well as the parton production during this period. 
It is well-known that there is a strong magnetic field from the colliding relativistic heavy ions. For 
the relevant overlapping region, it is mostly directed in the direction of the angular momentum 
which is perpendicular to the reaction plane as depicted in Fig. [TJ This has been studied in the 
Hadron-String-Dynamics (HSD) model [7]. For example, in Au + Au collision at y/S^N = 200 GeV, 
and impact parameter b = 10 fm, the magnetic field eBy at the center rises from zero at t —0.1 fm/c 
(N.B. t = is refereed to the time when the two nuclei are maximally overlapped.) to a maximum of 
~ 4.5m^ at t = 0.05 fm/c and quickly decays to ~ 0.7m^ at t = 0.15 fm/c [7J. On the other hand, 
due to the cancellation from the two nuclei, the electric field due to the Coulomb potential of the 
charges is weak in the overlapping region which we depict as Ec in Fig. [Tj For example, in Au + Au 
collision at ^/Snn = 200 GeV, b = 10 fm, the peak electric fields in the x-direction eE^ ~ 0.5 
at a: = 1 — 3 fm and t = 0.05 - 0.1 fm/c. This is an order of magnitude smaller than the peak 
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Figure 1: (color online) The reaction is taken to be along the x — z plane so that the angular 
momentum L and the magnetic field B in the center of the colliding nuclei are along the y axis. 
The electric field Ep (in red) induced by the transient magnetic field due to the Faraday law is 
circling around the y axis. The red arrows pointing outward indicate the electric field Ec due to the 
electric charges which is weak in the overlapping region because of the cancellation of the Coulomb 
potentials from the charges in the two colliding nuclei. 



magnetic field in the center of the overlapping region at t = 0.05 fm/c. However, we should point 
out that there is another transient electric field which has not been taken into account. Since the 
strong magnetic field decreases very quickly [3l [7j, i.e. eBy decreases maximally from ~ 4.5 at 
t = 0.05 fm/c to ~ 0.7 at t = 0.15 fm/c [7]. According to Faraday's law, this rapidly changing 
magnetic field leads to a transient electric field Ep which is circling the y-axis in the x — z plane 
(depicted in Fig. [T]) with a magnitude 

\eEF\ = = 19r(/m)m2, (4) 

where r in fermi is the distance from the y-axis. For h = 10 fm and the radius of the Au 
R = 1.2^4^/^ = 7 fm, the surface of the overlapping region around the reaction plane is at r = 3 
fm. At this distance from the center, {eEpl in the x — z plane is ~ 57 which is an order 
of magnitude larger than the maximum magnetic field and is even several times larger than the 
central NN force near the core at 0.7 fm of the Argonne V18 potential [8j. With such a strong, 
albeit short-lived, electric field Ep, the charged partons produced at this early stage of the collision 
(~ 0.05 fm/c) can acquire an impact with an initial momentum along the reaction plane as large 
as |Ap||| = \eEp\At = 560 MeV/c. This is comparable to the average p_L ~ 1 GeV/c [T3]. In order 
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to have a more realistic picture, one needs to fold in the timeline of the parton production during 
this initial phase. Using the Heavy Ion Jet INteraction Generator (HIJING) model to illustrate 
the space-time evolution of the parton production [T3], it is found that the partons are mostly pro- 
duced in the initial phase of the collision before 0.5 fm/c and the 'semi-hard' partons are produced 
between t = —0.15 fm/c and 0.15 fm/c. The total parton number produced as a function of time 
for the central collision of Au on Au at = 200 GeV is plotted in Fig. 3 in Ref. [13]. Before the 
nuclei overlap at t = 0, most of the partons are produced via initial state bremsstrahlung. Since 
these radiations are almost collinear, these partons move along the beam direction and have large 
rapidities [y is narrowly peaked at 2 at t = 0). These partons with large longitudinal momentum 
will leave the interaction region after semi-hard scattering and will not rescatter with the beam in 
leading twist. Since the STAR experiment of the azimuthal correlation are measured in the mid- 
pseudorapidity region {\ri\ < 1.0) and 0.15 < pt < 2 GeV/c, these partons produced before t = 
hardly will contribute to the experimental measurements. Between t = and t = 0.2 fm/c where 
the magnetic field eBy varies from 3 at t = to a maximum of 4.5 at t = 0.05 fm/c and 
down to 0.2 ml at t = 0.2 fm/c, the partons are produced via semi-hard scattering and the rate 
of production is the largest. The semi-hard scattering produces partons uniformly in the rapidity 
range y < 2 and fill up the mid and central rapidity region. These partons are copiously produced 
in the duration where the magnetic field varies the most and will be imparted large momentum 
from the induced electric field parallel to the reaction plane. Those produced at t ~ and afterward 
will see a net decrease of eBy and, thus, will experience a net electric force. Those on the surface 
of the overlap region, the maximum momentum transfer is 560 MeV/c as estimated above. After 
the semi-hard scattering, the final state radiation will also raise the mid and central rapidity of 
the produced partons, yet the magnetic field does not decrease much {eBy changes from 0.2 
at t = 0.2 fm/c down to near zero at t = 0.3 fm/c.) Thus, these partons will not be affected 
much by the electric force. But, the total partons produced after semi-hard scattering (i.e. after 
t = 0.2 fm/c) is small (~ 20% of the total) [13J. To sum up this part, the majority of the partons 
produced with rapidity y < 1 are produced via semi-hard scattering and it coincides with the large 
decrease of the magnetic field (from 3 at t = to 0.2 at t = 0.2 fm/c) where the induced 
electric field Ep gives a large momentum transfer in the direction parallel to the reaction plane. 
Moreover, the magnetic field, which is perpendicular to the reaction plane, will also exert a force 
on the moving charged partons along the reaction plane. In other words, both the electric field Ep 
and the magnetic field exert strong forces on the charged partons (quarks) along the reaction plane 
according to the equation F = qEp + qv x B. 



4 Charge-dependent Azimuthal Correlations 

We don't know how much the induced electric field will affect the subsequent hadronization and 
rescattering of the emergent charged hadrons. It will require dynamical modeling by incorporating 
the electromagnetic effect in this early stage as the initial condition for a more detailed description of 
the ensuing transport. In fact, there are calculations of the charge- dependent azimuthal correlations 
based on the Hadron String Dynamics (HSD) transport [7] and A Multi-Phase- Transport (AMPT) 
model P3] which found that the azimuthal correlations are much smaller than found experimentally. 
However, these calculations have not included the effect of the induced electric field from the 
Faraday law. It would be interesting to see if the induced electric field could have a large effect 
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in these calculations. Partly motivated by the presence of this large transient electric field which 
imparts momentum transfer to the partons along the reaction plane initially as discussed above, 
we shall make an ansatz that the azimuthal correlations mainly originate from the clusters on the 
surface of the overlapping region as discussed earlier and, furthermore, the positive- and negative- 
charged hadron pairs emerging outward from these surface clusters, as depicted in Fig. 2(a), will 
move mainly along the reaction plane. This parallel motion could be enhanced by the subsequent 
Coulomb electric field when the non-overlapping part of the charged nuclei fiy apart [7J. This is 
illustrated in Fig. 



2 a 



where the positive- and negative- charged hadron pairs emerging from a 
surface cluster in the positive x-hemisphere of the overlapping region with the positively charged 
particle momentum in the (x, z) = {+, ±, — ) direction and negatively charged particle momentum 
in the (x, y, z) = (+, ±, +) direction. But both have small azimuthal angles relative to the reaction 
plane with sin0 ~ e. On the other hand, the charged hadron pairs from a surface cluster which 
were originally produced from the partons heading toward the inner region will interact strongly 
with the denser medium [3] and emerge in the negative x-direction with large azimuthal angles as 



illustrated in Fig. 2(b) 
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Figure 2: (color online) (a) An illustration of the outward going charged hadron pairs which are 
mainly parallel to the reaction plane so that the azimuthal angles (j)^{out) and (f)^{out) are small, 
(b) Two charged pairs originating from one cluster with one pair going outward from the overlapping 
region and the other going inward through the overlapping nuclear medium. 

It is learned that the effective size of the cluster Mduster is ~ 3 in the Cu + Cu and Au 
+ Au collisions at center of mass energy of 200 GeV [lOj. The average number of pairs is 
{Npair{Npair — l)/2) = M^iusterf^ ^o that the total number of pairs is 



Nls + Ni 



us 



cluster M^l^g-^^^f Z, 



/2, 



(5) 



from the Poisson distribution. Here, Nls is the number of like-sign (LS) pairs and Nus is the 
number of unlike-sign (US) pairs. With Mduster ~ 3, the average sum of Nls and Nus is ~ 4.5 for 
each cluster. We shall consider the case from 4 up to 8 such pairs of Nls + Nus with the same 
number of positive- and negative- charged hadrons emitted from each cluster. 

Now let's consider the case of 2 positive and negative pairs from a cluster according to the 
ansatz - one pair in the positive x-direction with small azimuthal angles and the other one in 



the negative x-direction with large azimuthal angles in Fig. 2(b) For the like-sign (LS) pairs 



6 



(sin0-|-^i(out) sm0-|-^2(in))LS ~ 0(e) and (cos0±,i(out) cos0±,2(in))LS < in this case. These are 
in agreement with Eqs. ([T]) and (j2]). For the unhke-sign (US) pair, there are three terms for 
((+, l)(out), (-, l)(out)), ((+, 2)(in), (-, 2)(in)) and ((+, l)(out), (-, 2)(in)) and ((-, l)(out), (+, 2)(in)) 
pairs 



(sin 0+^i(out) sin(/)_^i(out))us 
(sin 0±,i(out) sin0±,2(in))us 
(sin 0+,2(in) sin0_,2(in))us 



0{e)- 



Sus > 0, 



where 



with /) 



us 



sin (f)a sin (pp fx 



us\ 



USWa, 



TT Jtt/2 

being the distribution of 0^ and 0^ for n/2 < 



(6) 



(7) 



< vr. Since sine is semi- 



positive definite, Sus is positive which is in agreement with Eq. ^ in sign. 
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Figure 3: (color onhne) Three pairs of positive- and negative-charged particles with one pair going 
outward and two pairs inward in (a) and two pairs going outward and one pair inward in (b). 



Similarly, 



(cos0+,i(out) cos0_,i(out))us ~ 1; 
(cos0-|-^i(out) cos0=p^2(in))us = 2c < 0; 
(cos0+,2(in) cos0_,2(in))us = Cus > 0, 



where 



2 

C = — d(f)a 
71 Jn/2 

Cus = (^) , 

Vr Jtt/2 



COS 0a(in)g(0«); 

cos 0Q (in) cos (in) fus ( 



(9) 
(10) 



Both angles 0a(in) and 0^(in) are between 7i/2 and vr; therefore, Cus > 0. To the extent that 
(1 -|- 2c -|- Cus) > 0, it would be in agreement with Eq. ^ in sign. Taking a crude estimate that 
the inward particles emerging in the negative x-hemisphere have a uniform distribution and 0o(in) 
and 0^(in) are uncorrelated, then c = —2/7i and Cus = 4/(7r)^. In this case, (1 + 2c + Cus) = 
(1 -2/7r)2 > 0. 
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We can extend this scenario to more general cases with equal numbers of positively /negatively 
charged particles (M_|_ = M_ = M) in the positive x-hemisphere and similarly with A^_|. = = 
charged particles in the negative x-hemisphere, all originate from one cluster. For the LS case 

{sin (pa sin (p ,3) Ls oc M{M - l)e'^ + 2MNe s + N{N - 1)Sls, (H) 

where s = J^^2 sin0Q,(in)g(a) gives the contribution from the cross term with one particle in the 
-\-x hemisphere with sin 0(o-ut) ~ e and the other one in the —x hemisphere with a distribution of 
g{a). Sls is the same as Sus in Eq. ^ except the distribution function fus is replaced with fus- 
Similarly, 

(cos 0,, cos 0/3) L5 OC M {M - 1) + 2M N c + N {N - 1)Cls , (12) 



where Cls is the same as Cus in Eq- (10) except the distribution function fus is replaced with fus- 
On the other hand, for the unlike sign case 

{sin (pa sin (f)p) us oc M^e^ + 2MNe s + N^Sus, (13) 
{cos (pa cos (pp) us oc + 2MNC + N^Cus- (14) 

We give a table of these results for the cases with a total number of charged particles K {K = 
2(M + N)) being from 4 to 8 in each cluster in Table [l] 

Table 1: Contributions to charge-dependent azimuthal correlations for 4 to 8 charged particles from 
2 to 4 pairs of positive- and negative-charged hadrons. 



Total Numbers 


M 


N 


(sin (pa sin (pp)Ls 


{cos (pa cos (Pp)ls 


(sin (pa sin (pp)us 


{cos (pa COS (pf})us 


(Charged Particles) 














4 


1 


1 


0{e) 


2 c 


Sus 


1 + 2c + Cus 


6 


1 


2 


Sls 


Ac + 2 Cls 


A Sus 


1 + 4c + ACus 




2 


1 


0{e) 


2 + 4c 


Sus 


4 + 4c + Cus 




Sum 


6 


Sls 


2 + 8c + 2 


5 Sus 


5 + 8c + 5Cus 


8 


1 


3 


6 Sls 


6c + QCls 


'9 Sus 


1 +6c + 9Cus 




3 


1 


0{e) 


6 + 6c 


Sus 


9 +6c + 9Cus 




2 


2 


2 Sls 


2 + 8c + 2 


A Sus 


4 + 8c + 4 Cus 




Sum 8 


8 Sls 


8 + 20c + 8 Cls 


14 Sus 


14 + 20c + 14 Cus 



Due to the reflection symmetry with respect to the x-axis, the charged hadron pairs from the 
—X part of the overlapping region have the same azimuthal correlations. The total azimuthal 
correlations for the like-sign charged particles are 

(sin(0,)sin(0;3))L5 = ^[P{A-X)0{e) + P{6-X)SLs + 8P{8,X)SLs], (15) 

(cos(0„) cosi^p)) LS = |[P(4, A) 2c + P(6, A) (2 + 8c + 2Cls) + ^(8, A) (8 + 20c + 8Cls)], 

(16) 
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where Z is the normahzation factor and PiK-, A) is the Poisson distribution for K particles with a 
mean of A. 

Likewise, the total azimuthal correlations for the unlike-sign charged particles are 
(sin(0,)sin(0;3))c/5 = ^[P(4,A)5t;s + 5P(6,A)5^5 + 14P(8,A)5f;5], 

(cos(0„) cos(0^))[;5 = ^ [P(4, A) (1 + 2c + Cus) + P(6, A) (5 + 8c + W^s) + P(8, A) (14 + 20c + 14C 



There are few prominent features that we can readily observe from Eqs. ( 15 )-( 18 ) . With the mean 
value being the average cluster size A = Mduster = 3, the Poisson probabilities are P(4; 3) = 0.168, 
P(6; 3) = 0.0504, P(8; 3) = 0.00810 (We have ignored the 10-particle contribution due to its small 
P(10; 3) = 0.00081 which is an order of magnitude smaller than P(8;3)). The correlations are 
dominated by the case with 4 particles. The noted features are: 

• Since the 4-particle contribution to (sin(0o) sin(0/3))i5 is negligible (i.e. 0(e)), it is expected 
to be small as compared to other correlations which is consistent with Eq. ([T|. 



Since P(4; 3) is the leading term and the fact that c < 0, (cos(0a) cos(0/3))l5' in Eq- (16) turns 
out to be negative. This is consistent with experiments in Eq. ([2]). 



• We see from Eqs. (18) that if c ~ —0.5 which is close to the uncorrelated estimate of — 2/7r, 
then the out-out correlations are largely canceled out by the out-in cross correlations for each 
of the 4-, 6-, and 8-particle cases. Furthermore, we notice that the coefficients of Sus and 
Cus are identical in each case with different particle numbers. Thus, if Sus ~ Cus which is 
true for the flat distribution of fus in Eq- ([T)), then the sine and cosine correlations of the 
unlike-sign pairs are about equal, i.e. (sin(0a) sin(0^))[/5 ~ (cos(0q) cos(0^))[/5 which is what 
is found experimentally in Eq. (|3|. 

Thus, the present model with simplistic scenarios and assumptions seem to be able to capture the 
generic features of the azimuthal correlations in Eqs. ([T]) - ([3]). 

To be more specific, let's take c = —0.55 and Cls = Sus = Cus = 0.8 and A = Mduster = 3, we 
find 

(cos(0c.) cos(0/3))l5 : (sin(0o) sm{(f)p))us ■ (cos(0o) cos{(l) ,3)) u s = -0.46 : 1 : 1.09 (19) 

These ratios are in quantitative agreement with those of experiments in the centrality range of 35% - 
65% as illustrated in Fig. 2 in Ref. [Bj. The fact that (sin(0cj) sin(0^))[75 and (cos(0q,) cos(0/3))c/5 are 
near identical is understood as due to the cancellation between the out-out and out-in correlations 
for the unlike-sign charged particles. The fact that (sin((^a) sin(0/3))c/5 is negative is due to the 
fact that it is dominated by the out-in cross contribution (i.e. terms associated with c) which is 
negative. 



As for the smallness of (sin(0Q,) sin(0/3))i,5', we see that the ratio between Eqs. (15) and (17) 
depends only on |^ 

(sin(0jsin(0^))z.5 ^Q o^g/'^Lg. ^20) 
(sin(0c.)sin(0^))c/5 ' Sus ' 
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If is about 0.5, it would roughly explain the experimental finding that ^sin(^") sin(^^))^g about 
an order of magnitude smaller than (sin(0a) sin(0^))[/5 and (cos(0a) cos(0/3)){/5 at 60% centrality. 

5 Conclusion 

To conclude, we find the cluster model description with the assumption that the charge-dependent 
azimuthal correlations originate mainly from the clusters on the surface of the overlapping region 
of the colliding nuclei and that the outward going charged hadron pairs from these surface clus- 
ters are parallel to the reaction plane can describe the azimuthal correlations of both the like-sign 
charged pairs and unlike-sign charged pairs in a qualitative manner. In particular, we find that 
the hard to understand near equality of the sine and cosine correlations of the unlike-sign pairs 
i.e. {sm{(j)a) sm{(j)j3))us — (cos(0q,) cos(0^));75 is due to the cancellation between the out-out pairs 
and the out-in cross correlated pairs. In addition, the centrality dependence of the correlations is 
naturally explained in terms of the cluster model [lU] and the surface cluster dominance is suggested 
from the suppression of the cluster size of the near-side pairs as compared to the away-side pairs 
in experiment |lOj. We also point out that there is a strong transient electric fields induced by the 
rapidly varying magnetic field in the initial phase of the collision. Folding in the timeline of parton 
production, we find that most partons produced on the surface during the semi-hard scattering after 
t = with mid rapidities coincide with the rapid decrease of the magnetic field and will experience 
a momentum transfer as large as several hundred MeV from the induced electric field along the 
reaction plane. This could have some relevance to our assumption of the outward going charged 
hadron pairs. To verify this simplified picture presented here and explain the data more quanti- 
tatively, one needs to incorporate the large transient magnetic field and its collateral electric field 
from the Faraday law in the very early stage of parton production (0.1 fm/c) in dynamical models 
such as Hadron String Dynamics (HSD) [15j, HIJING \16\, AMPT [I7|, NeXus yLS], momentum 
kick model [19], and 3D Relativistic Hydrodynamics [20]. Only when the background contribution 
from the parity-non-violating QCD + electromagnetic effects are understood to enough precision 
will one be able to address the existence of the chiral magnetic effect in these azimuthal correlations. 

This work is partially support by U.S. DOE Grants No. DE-FG05-84ER40154. The author 
wishes to thank H. Huang, F. Wang, C.Y. Wong, and N. Xu for useful discussions and particularly 
C.Y. Wong for reading the manuscript. He also thanks Mingyang Sun for drawing the figures for 
this manuscript. 

References 

[1] B.I. Abelev et al. [STAR Collaboration], Phys. Rev. Lett. 103, 251601 (2009); B.I. Abelev et 
al. [STAR Collaboration], | arXiv:0909.1717] . 

[2] D. Kharzeev, R.D. Pisarski, and M.H.G. Tytgat, Phys. Rev. Lett. 81, 512 (1998); D. Kharzeev, 
Phys. Lett. B 633, 260 (2006); K. Fukushima, D.E. Kharzeev, and H.J. Warringa, Phys. Rev. 
D78, 074033 (2008). 

[3] D.E. Kharzeev, L.D. McLerran, and H.J. Warringa, Nucl. Phys. A803,227 (2008). 



10 



[4] LB. Khriplovich and A.S. Rudenko, |arXiv:1003.3339|. 

[5] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 92, 062301 (2004); N. Borgliini, 
RM. Dinh, and J.Y. Ollitrault, Phys. Rev. C 64, 054901 (2001); A.M. Poskanzer and S.A. 
Voloshin, Phys. Rev. C 58, 1671 (1998). 

[6] A. Bzdak, V. Koch, and J. Liao, Phys. Rev. C81, 031901 (2010). 

[7] V. Voronyuk, V.D Toneev, W. Gassing, E.L. Bratkovskaya, V.P. Konchakovski, and S.A. 
Voloshin, |a rXiv:1103.4239v 2 [nucl-th]]. 

[8] R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, Phys. Rev. C 51, 38 (1995). 

[9] K. Eggert, et. al., Nucl. Phys. B86, 201 (1975); B. Alver et al. (PHOBOS GoUaboration) , 
Phys. Rev. C 75, 054913 (2007). 

[10] B. Alver et al. (PHOBOS GoUaboration), Phys. Rev. C 81, 024904 (2010). 

[11] B.B. Back et al. (PHOBOS GoUaboration), Phys. Rev. C 74, 011901 (R) (2006); B. Alver et 
al. (PHOBOS GoUaboration), Phys. Rev. C 75, 054913 (2007). 

[12] F. Wang, Phys. Rev. C81, 064902 (2010). 

[13] K.J. Eskola and X.N. Wang, Phys. Rev. D49, 1284. 

[14] G.L. Ma and B. Zhang, Phys. Lett. B700, 39 (2011). 

[15] W. Ehehalt and W. Gassing, Nucl. Phys. A602, 449 (1996); W. Gassing and E.L. 
Bratkovskaya, Phys. Rep. 308, 65 (1999). 

[16] X.N. Wang and M. Gyulassy, Phys. Rev. D44, 3501 (1991). 

[17] Z.W. Lin, G.M. Ko, B.A. Li, B. Zhang, and S. Pal, Phys. Rev. C72, 064901 (2005). 

[18] H.J. Drescher, P.M. Liu, S. Ostapchenko, T. Pierog and K. Werner, Phys. Rev. C65, 054902 
(2002). 

[19] G.Y. Wong, Phys. Rev. C84, 024901 (2011). 
[20] Y.Hama, et al., Nucl. Phys. A774, 169 (2006). 



11 



